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Abstract 

The main objective of this work is the pursuit of a detailed analysis of the cosmological solutions 
corresponding to each fixed point belonging to an exclusive form of f(R) initially introduced by Hu 
and Sawicki in the background of Palatini f(R) gravity. Fixed points and their associated solutions for 
the anisotropic Bianchi type I and V models are derived by utilizing an approach used to solve 
autonomous ordinary differential equations. Building on this background, an exhaustive study on the 
phase space is accomplished and the existence of a sequence of evolution that starts from unstable 
anisotropic universe to a stable phase of late times de Sitter expansion through the standard radiation 
and matter dominated era has been established. 


1. Introduction: 

In recent years, modified theories of gravity have earned great interest among researchers all over the globe as these 
theories are capable of reproducing the accelerated phase of the late universe. Although many theories of modified 
gravity have been constructed to explain such late time acceleration,/^) gravity [1-10] is considered one of the most 
extensively studied and as the simplest of such theories. In f(R) theories, the Ricci Scalar R in the Lagrangian density 
corresponding to the standard Einstein-Hilbert action is replaced by some arbitrary function of R. The field equations 
corresponding to the modified action can be obtained executing different formalisms of f(R) gravity and in this analysis 
we pay our attention to the Palatini formalism. In this formalism, field equations are obtained by varying the action 
corresponding to the metric as well as the connection, considering them as two independent variables. An additional 
assumption in Palatini formalism is that the matter action is independent of the connection. In addition to the fact that 
the field equations correspond to mathematically simple second order differential equations, the Palatini formalism 
has gained sufficient interest from the perspective of the cosmological tests [11]. Palatini formalism of f(R) gravity is 
capable of reproducing the three post inflationary epochs of the standard cosmology for some of the cosmological 
models based on a power law functional form, which is not possible in case of the standard metric formalism. In this 
analysis, we study an exclusive form of f(R\ as originally proposed by Hu and Sawicki [11-18] within the Palatini 
approach. The cosmological consequences of Palatini/fR) gravity for an exclusive form of f(R) initially introduced by 
Hu and Sawicki in the background of Friedmann-Leimaitre-Robertson-Walker (FLRW) model has been investigated 
in [11]. It has been observed that the Hu-Sawicki model is indistinguishable from the standard ACDM scenario and is 
able to accomplish the sequence of radiation dominated, matter dominated and accelerated expanding period in the 
Palatini formalism. Hu and Sawicki have provided the metric variation for a class of f(R) gravity which accelerates the 
expansion without considering a cosmological constant and also satisfies the cosmological as well as solar system tests 
corresponding to the small field limit of the parameter space [12]. A detailed analysis of the dynamics corresponding to 
Hu-Sawicki f(R) model is presented in [13]. A full phase space analysis has been performed for the FFRW metric and 
several de Sitter fixed points have beenidentified along with a matter like unstable fixed point for a particular value of 
the model parameters. Authors of [14] have focused on the analysis of the viable Hu-Sawicki f(R) model, which can 
reproduce the late time accelerated expansion era and also can recover the result that is consistent with the General 
Relativity (GR) on local scales. A combined investigation on the theoretical as well as observational constraints which 
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impart some important information about the viable parameter space of this explicit form of f(R) model have been 
provided. Authors of [15] have explored the two viable/(R) models (Hu-Sawicki and Nojiri-Odinstov), by 
investigating their corresponding cosmological evolution and studying their scalar-tensor representation. They study 
the scalar-tensor equivalence of both the viable/(R) models, which provides a very useful tool to extract the necessary 
information about the possibility of occurrence of cosmological singularities. 

Despite being mathematically simpler in comparison to the other formalisms of f(R) gravity, the system of 
field equations corresponding to Palatini /(R) gravity are non linear in nature and to study the behaviour of the 
non linear system we use a suitable technique known as Dynamical System Approach (DSA) [19-23]. Such a 
technique enables us to study the behaviour of the equilibrium points corresponding to the equations of the 
dynamical system. Fay et al have studied the late time acceleration for FLRW model in Palatini /(R) gravity 
adopting DSA [5]. 

The acceleration phase of the late universe can be studied in the background of Bianchi type I and V models 
which are the simple generalization of the FLRW model. FLRW model is spatially homogeneous and isotropic 
[5, 24-26] whereas Bianchi type I and V models are spatially homogeneous but anisotropic [27-30]. The 
simplest representation of isotropic and spatially homogeneous stage of present universe can be excellently 
analysed in the background of FLRW model. However, the early stage of the Universe, which could have been 
anisotropic [3 1-37], cannot be explained in the background of this isotropic model. It maybe noted that the 
Cosmic Background Explorer (COBE) confirmed the existence of anisotropy in 1992 [31] and for this discovery 
the Nobel Prize in Physics in 2006 is awarded to J C Marther and G F Smoot. In addition, some recent probes and 
experiments, for example, Wilkinson Microwave Anisotropy Probe (WMAP) [32- ] and Planks results [36,37] 
also have confirmed the presence of small anisotropies. In a recent paper [38], the authors have used two 
methods, namely the hemisphere comparison (HC) method and dipole fitting (DF) method in probing the 
possible preferred direction in the distribution of type la supernovae (SNIa). They used these methods in order 
to test cosmic anisotropy with supernovae data. Evidences of cosmic anisotropy are still accumulating all over 
the globe. The existence of cosmic anisotropy cannot be fully investigated in the background of isotropic FLRW 
cosmology. This provides a strong motivation to look for some alternative models to explain such anisotropy. 
Bianchi models [39-53] pave a way for the researchers to describe the anisotropic phase of early universe and the 
acceleration phase of late universe. 

Our analysis has been carried out for an exclusive form of/(R), as originally proposed by Hu and Sawicki, 
within the Palatini approach for anisotropic Bianchi I and V models. The main motivation to use the Palatini 
formalism in the Hu-Sawicki scenario is to check whether such a model is able to reproduce the sequence of 
cosmological evolution in the presence of anisotropic models or not. Although the FLRW cosmology within the 
Palatini approach [11] and the anisotropic Bianchi III cosmology within the metric approach [52] in the 
background of Hu-Sawicki /(R) model has been able show the sequence of the evolution of the Universe, the 
Palatini version of Hu-Sawicki /(R) gravity has not been yet investigated in presence of anisotropic models. Such 
anisotropic models are considered as a preferred area of research as it can explain the anisotropy prevailing in the 
early stage of the Universe. 

The content of the present work is outlined as follows. Section 2 presents the field equations of f(R) gravity in 
the Palatini variational technique. In section 3, we discuss the Hu-Sawicki /(R) model and its cosmological 
viability. In sections 4 and 5, we derive the fixed points along with their solutions, stabilities and also pursue the 
phase space analysis for anisotropic Bianchi I and V cosmologies, respectively. Ultimately, section 6 contains the 
conclusion of the work. 


2. The Palatini approach in f(R) gravity 


The/(R) gravity action that is defined as the modification of Einstein-Hilbert action is written as 


5 = 


fd*xj=g 


~~yf W + L m + L r 
,2k l 


( 1 ) 


where/(R) denotes a function of the Ricci scalar R, g represents the determinant of metric tensor L r and L m 

are defined as the Lagrangians of matter and radiation respectively, n 2 = 87 rG and G is the gravitational 
constant. As disscussed in [5], matter and radiation are introduced ad-hoc into the equations of motion. 

In Palatini approach, the basic idea is to consider the metric g^ v as well as the affine connection T^ w as two 
independent variables. The field equations arising due to the variation of the action (7) corresponding to the two 
independent variables are 

m v - J/(% F = (2) 


2 
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Va (j=gfmgn = o 


(3) 


where the prime denotes the derivatives with reference to R and Vdenotes the covariant derivative. T^ is the 
tensor of energy momentum defined as 


2 Hlr L r ) 


(4) 


3. The Hu-Sawicki/(i?) model 


In this analysis, particular attention is given to a unique form of f(R) model known as Hu-Sawicki f(R) model 
[11-15,52], given by 


f(R) = R 






(5) 


where c ly c 2 , n are constant parameters and the mass scale is taken as m 2 = K 2 p/3. A detailed investigation of the 
cosmological dynamics corresponding ton = 1, in case of Hu-Sawicki f(R) model are presented in [1 1,13, 52]. 
Following [ 11 , 13, 52] , we perform our analysis considering n = 1, as the value of n is totally unconstrained by 
the present cosmological data. Thus, for n = 1, the above equation reduces to 

2A 

f(R) = R- - - (6) 

Y + 1 

where A = m 2 q / 2c 2 and p 2 = m 2 /c 2 are two dimensionless parameter, while m 2 = Hq D m0 ( H 0 and D m0 are 
the present-day values of Hubble and matter density parameters, respectively). In the regime R <C p 2 the 
correction to GR is found to be negligible [14]. The above model is found to be practically indistinguishable from 
the ACDM model in the regime R [i 2 [1 1-14]. Authors of [12] suggest that R /i 2 , even in the present 

epoch. In our analysis, we have considered the condition R /i 2 . 

A ACDM model is the simplest cosmological model which can drive the Universe into the accelerating de 
Sitter phase. The cosmological constant A was originally introduced by Einstein in GR to keep the Universe 
static. But, now a small positive value of the cosmological constant is considered as a possible candidate for dark 
energy. The ACDM model has provided some well defined predictions that have been able to withstand the 
rapid as well as continuous advancement of cosmological observational tests. Observations of the large scale 
structure and recent accelerated expanding universe are in excellent agreement with the ACDM model. Despite 
these major advantages, ACDM model has several problems. One of such problems is the cosmological constant 
problem, which has not been satisfactorily resolved until today. [1 1, 12, 14] have shown that Hu-Sawicki f(R) 
model can explain the accelerated expansion without considering a cosmological constant. 

Let us now use the scalar-tensor theory of f(R) gravity to examine the sudden singularity for the Hu-Sawicki 
model. 

The action corresponding to the scalar-tensor theory of f(R) gravity is written as 


S = 


fd 4 xj=g 


—(0R - Vm + L m + L r 
2k 2 


(7) 


where the scalar field </> and its potential V(<p) is given by 

</>=/', V(<fi) = Rf — f 


( 8 ) 


In general the explicit expression of the scalar potential in terms of the scalar field for a general n for Hu- 
Sawicki/fR) model is not possible to achieve. However, this is possible corresponding to the case n = 1, such 
that the scalar potential is given by [14,15] 

V(</>) = m 2 (1 ~4>)± 2V2A^ 2 (1 - 0) + 2A (9) 

Substituting the value of /./, 2 and A in the above equation we get 


= H 2 n 2 n 


ci + (1 - <j>) ± 2^ci(l - <j>) 


( 10 ) 


It is observed that the sudden singularity, where R —> oo, exists for 

0-i, v^H 2 n 2 m0 ^ (ii) 

C2 


3 
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Figure 1. Evolution of the scalar potential (V±(0)) corresponding to Hu-Sawicki model for n = l,Ci = 2 ,c 2 = 0.1andi2 m0 = 0.27. 


The evolution of scalar potential V((f)) for n m 1 is plotted in figure 1 . This figure shows that the two branches 
of the scalar potential ( V± (</>)) consist of different asymptotically stable points. While the upper branch (V + (</>)) 
stops at the singular point 0=1, the lower branch (VL(0)) accomplished the asymptotically stable de Sitter 
evolution (see [14]). Therefore in order to avoid the singularity, we have to stay on the lower branch of the scalar 
potential ( V-(<j))) which leads to the condition 

%=D < H 0 2 ^A (12) 

Thus, the singularity can be avoided depending upon the initial conditions and the model parameter values. 


4. Bianchi I metric in Palatini/CR) gravity 


Our analysis begins with one of the simplest type of anisotropic models known as Bianchi type I. The line 
element of this anisotropic model is defined as 

ds 2 = —dt 2j r A 2 (t)dx 2 + B 2 (t)[dy 2 + dz 2 ] (13) 


where A(t) and B(t) are considered as the expansion scale factors. x,y, z represent the comoving coordinates 
whereas t denotes the cosmological time. 

The generalised Friedmann field equation in case of the anisotropic Bianchi type I model can be written as 


V 

3 


e + 


it 

2 r 


- 2fcr 2 -/= n(p m + 2p r ) 


(14) 


where a dot represents differentiation with respect to cosmological time. 0 and a stands for the volume 
expansion scalar and the shear scalar, respectively and are defined as 


. d A 6 

6 = 3— = -h 2— 

a A B 


(15) 


3 [a 


12 


(16) 


where a acts as the average scale factor. p m and p r appear as the matter energy density and radiation energy 
density, respectively, which satisfy the conservation equations given by 


Pm + & Pm ~ 0 
Pr + ~®Pr = 0 


4 


(17) 

(18) 
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The trace free Gauss-Codazzi equation in case of Bianchi type I model for Palatini/fR) gravity is defined as 


& = 


f"R 

f 


a 


(19) 


The above equation is very beneficial in describing the evolution of shear with reference to time. Considering 
the trace of equation (2) and using the conservation equation satisfied by the energy density of matter (17) in the 
resulting trace equation, one obtains 


R = 


2 Opm 
f"R - f 


= -e 


f'R - 2 f 


( 20 ) 


Using equation (20) in the generalised Friedmann field equation (14), one can obtain 
a2 __ 6f r <J 2 + 6k 2 (p m + p r ) + 3 {f'R — f) 

0 -W- 

where 

, , . 3f'R ] S fiTR-Tf) 

? 2 f e 2 f(f'R — 2/') 


( 21 ) 


( 22 ) 


To perform the stability as well as the phase space analysis, it is essential to construct a system of first order 
differential equations of the cosmological equations using the DSA technique. In this approach, the analysis 
begins with the introduction of appropriate dimensionless variables. We define the following suitable variables: 


V3cr 


tt r = 


3 ifPr 

fee 2 


fi» 


ef£m r 

fee 1 ’ 


3 (f'R — /) 
2 fee 2 


(23) 


The above variables are constrained by the modified Friedmann field equation (21) as 

1 = Yj 21 T U r T ~F x 


(24) 


Using the variables (23) along with the constraint equation (24), one can write the following evolution 
equations corresponding to the variables E, U r and % as: 

V' 

— = —[-3 + 3£ 2 +Q r -3x- 9C(R)x + 9D(R)(S 2 - 1)] 
dr 2 

— = 12 r [-l + 3£ 2 + fi r - 3x - 9C(R)x + 9D(fl)£ 2 ] 
dr 

dx 

— = x[3 + 3E 2 + f2 r - 3* + 9CfR)(l - x) + 9 D(R)Z 2 ] (25) 

dr 


where r = In a 

The expression for C(R) and D(R) are given as 


1 Rp 


3d (f'R - f ) 


D(R) 


1 f'R 
30 — 


1 ( f'R - 2f) f'R 

1 (fR - 2f) f" 

3 (f'R-f)f 


(26) 

(27) 


In the dynamical system (25), it is realized that the only time (or r) dependence is contained in the 
parameters C{R) and D(R). In our analysis, we have considered the condition R p 2 as the Hu-Sawicki model 

is found to be practically indistinguishable in this regime [11-14]. The condition R <C p 2 is avoided as the 
correction to GR in this regime is found to be negligible [14]. In the regime R p 2 corresponding to the 

Hu-Sawicki model, one can find that 


C{R) 


2/i 2 

3 R 


D(R) = -|q 
3 



(28) 

(29) 


Therefore in the regime R p 2 ,C(R) landD(R) <C 1 and these conditions are in general valid for all 
the heteroclinic orbits in the dynamical system. 

As C(R) < 1 and D(R) <C 1, they can be neglected in the evolution equations (25) as compared to the other 
terms As a result the dynamical system (25) does not depend on time (or r) explicitly, which proves that the 
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Figure 2. Evolution of EoS and deceleration parameter along with the dimensionless dynamical variables for Bianchi I cosmology in 
the background of Hu-Sawicki model. 


system we have obtained in case of Bianchi I cosmology is an autonomous system. Authors of [54] have provided 
a detailed investigation on the autonomous dynamical system in case of f(R) gravity. 


4.1. The fixed points and exact solutions: 

It is important to find the fixed points to study the cosmological behaviour of the dynamical system. The fixed 
points which can be used to procure the exact solution for the anisotropic Bianchi models are obtained by setting 
equation (25) equal to zero. In case of anisotropic Bianchi type I model, we encounter anisotropic as well as 
isotropic fixed points, given by 

P*i (E, fir, x) = (±1, 0, 0) 

P r : (E, fi r , x) = (0, 1, 0) 

Pm- (S, fi r , *) = (0, 0, 0) 

Pd- (E, fi r , X) = (0, 0, 1) 


Next, we proceed to describe the equation of state (EoS) parameter w e ff and the deceleration parameter q in 
terms of the dynamical variables as 


W eff = 


Peff E 2 + ifi r - (1 + 3 C(R))x + 3D(.R)(1 + E 2 ) - E 2 £ 2 + 2^ 


Peff 


q — — 1 + — 
2 


1 - £ 2 £ 2 

1 + E 2 + _ (i + 3 C(R))x + 3D(K)(E 2 + 1) + 2-4 


(30) 

(31) 


Inserting the values of each variable in equations (30) and (31) and considering C(R ) <C 1, D(R) <C 1 and 
£ « 1, one can find the EoS parameter w e ff and deceleration parameter q for each fixed point. The evolution of 
EoS parameter and deceleration parameter along with the dimensionless dynamical variables are plotted in 
figure 2. These figures show that the Universe starts from the anisotropic state (w e jy= 1/3, q = 2) that 
isotropizes with the standard radiation (w e jy= 1/3, g = 1) and matter dominated phase (w e jy = 0, q = 1/2) and 
finally stabilises to the late times de Sitter expansion epoch (w e jy = —l,q = —1). 

The relation between the deceleration parameter q , and the volume expansion scalar 0, is written as 

¥ = 4 ° + q> (32) 


Integrating equation (32) for the fixed points corresponding to q ^ — 1 and setting the big bang time t 0 = 0, 
we obtain 


e = —-— = ^ 

(1 + q)t t 


(33) 


We again integrate the above equation in order to find the solution of the average scale factor and the 
solution is found to be 
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Figure 3. Evolution of scale factor (a) corresponding to Bianchi I cosmology in the background of Hu-Sawicki model for a 0 = 0 O — 1. 


Table 1. Fixed points and their solutions for EoS parameter, deceleration parameter, average scale factor, shear and 
physical behaviour for Bianchi I cosmology in the background of Hu-Sawicki model. 


Points 


EoS 

parameter 

(W eff ) 

Deceleration 

parameter 

(q) 

Average 

scale 

factor(n) 

Shear 

O) 

Physical behaviour 

Pf 

0 

1 

3 

2 

a Q \t\\ 

cr 0 flo“ 3 |t| _1 

Decelerated expansion 

P r 

0 

1 

3 

1 

a 0 \t\i 

0 

Decelerated expansion 

Pm 

1 

0 

1 

2 

a 0 \t\i 

0 

Decelerated expansion 

Pd 

0 

-1 

-1 

a 0 \t\> 

0 

Accelerated expansion 


a = a 0 |t| c (34) 

where a 0 is a constant and £ = (1 + q)~ x 

For the fixed points having q = — 1, it is found that 0 = 0. Solutions of such fixed points corresponds to de 
Sitter solutions which are written as 

a = a 0 e^ eot (35) 

Therefore, we find that power law solution exists for q ^ — 1, whereas de Sitter solution exists for q = — 1 
and these solutions are plotted in figure 3. The average scale factor a, along with the EoS parameter w e ff, 
deceleration parameter q and physical behaviour for the fixed points are summarised in table 1 . 


4.2. Evolution of shear: 

In the following, we study the evolution of shear for the anisotropic fixed points. Such an evolution can be 
studied from the trace free Gauss-Codazzi equation (19). Equation (19) can be rewritten in terms of dynamical 
variables (23) as follows 


a = 


-= [1 + 3 D(R)^]0 2 

V3 


(36) 


For the anisotropioc fixed points P*, integration of the above equation leads to 

a = tT 0 fl _3 = cro^crkr 1 ( 37 ) 

The above equation indicates that the shear for the anisotropic fixed points evolves inversely with time and it 
is shown in figure 4. 


4.3. Stability analysis of the fixed points: 

In this section, we determine the eigenvalues of the system of equation (25), which are essential in studying the 
stability. The eigenvalues obtained by linearising the corresponding system of equations are as follows: 
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Figure 4. Evolution of shear (<r) corresponding to Bianchi I cosmology in the background of Hu-Sawicki model for a 0 = 1. 


Table 2. Eigenvalues and stabilities of the fixed 
points for Bianchi I cosmology in the background of 
Hu- Sawicki model. 


Points 

Eigenvalues [A b A 2 , A 3 ] 

Stability 

Pf 

[2,3,6] 

Unstable 

Pr 

[-1,1,4] 

Saddle 

Pm 

[-f, -1, 3] 

Saddle 

Pd 

[-4, -3, -3] 

Stable 


Py [A b A 2 , A 3 ] = [2 + 9I)(R), 3 + 9D(R), 6 + 9C(R) + 9D(R)\ 

P r : [A b A 2 , A 3 ] = [l, 4 + 9 C(R), -1 - jD(R ) 

Pm- [Ai, A 2 , A 3 ] = f-1, 3 + 9C(R), -- - —D(R) 

L 2 2 J 

Pd'. [A b A 2) A 3 ] = f-3 - 9 C(R), -4 - 9C(R), -3 - -C(R) - -9D(R) 
l 2 2 

In the regime R fi 2 corresponding to the Hu-Sawicki model, it is found that C(R) <C 1 and D(R) <C 1. 
Therefore, the eigenvalues of the Bianchi I model reduces to the form as listed in table 2. In the corresponding 
table we also summarise the stability of the fixed points of the dynamical system. For the anisotropic fixed points 
P±, we have all the three eigenvalues positive. Therefore, they correspond to unstable nodes. The points 
representing the radiation and matter dominated era given by P r and P m , respectively, are found to be saddle, as 
both have positive as well as negative eigenvalues. Finally, for the fixed point all the eigenvalues are found to 

be negative, which describe the stable de sitter expansion of the late universe. 


4.4. Phase space analysis: 

Let us now further investigate the system using phase space analysis in order to study the dynamical behaviour of 
the Bianchi I cosmology corresponding to the Hu-Sawicki model. The phase portrait analyses of the system are 
shown in figure 5. We can realize heteroclinic trajectories of the form 

Pf - Pr - Pm -> Pd 

This figure displays the progression of the initial unstable anisotropic phase to stable de Sitter expanding 
universe following the standard radiation and matter era. 
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5. Bianchi V metric in Palatini / ( R ) gravity: 


In the following, we study the more complicated form of anisotropic models known as Bianchi type V. The line 
element of this anisotropic model is defined as 


ds 2 = -dt 2 + A 2 (t)dx 2 + e~ 2 P x [B 2 (t)dy 2 + C 2 (i)dz 2 ] 


(38) 


where A(t), B(t ) and C(t) are considered as the expansion scale factors, x, y, z represent the comoving coordinates 
whereas t denotes the cosmological time andp is a constant. 

For this anisotropic model, the field equation (2) modifies to the following forms: 

r '\ 2 

+ 


/' 

/' 

/' 

f 


3/' 3 (/' 


2/' 2 If 


A 

A 


'MA 

2 f V A 

l£(A 

2fU 

if if (A 

2 f 2 f'\A 


if 
2 f 

if 
2 f 


+ 3— 
A B 


A B 
+ — 
A B 


(A 

B C) 

- + 

-1- 

u 

B CJ 

c) 

1 AB 

+ — 

H-F 

C) 

1 AB 

c ) 

| AB 

+ — 

+ - + 

Cj 

1 AB 

A) 

| AC 

\- 3- 

+-+ 

c) 

1 AC 

■\ 2 A 

B C" 

1 A 

B C 


1 


AC 

AC 

BC 

~BC 


A 

A 

B 

B 


2p2 

A 2 

2P 2 ' 

A 2 


BC C _ 2P 2 
BC + C A 2 


p = 0 


- /= 


-/ = kP 

2 


-■/= 

2 


- -/= kP 
2 


Equation (43) leads to the following equation 


^ A B C 
2 - = - + - 
ABC 


(39) 

(40) 

(41) 

(42) 

(43) 

(44) 


Solution of the above equation takes the form 

A 2 = cBC (45) 

where c is an integration constant. 

The generalised Friedmann field equation in case of the anisotropic Bianchi type V model can be written as 

h l 2 

- 2fa 2 - 6 p 2 f'( 2 K) -f=K(p m + 2p r ) (46) 


2 f 


e + lL 


2 f 


where 0, <j and 2 K stands for the volume expansion scalar, shear scalar and Gauss curvature of the 3-spheres, 
respectively, are defined as 


« = 3 i = 4 


B C 

— + — 
B C 


= 3 l = l(* + t) 

A 2\B C) 


(47) 
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(48) 

(49) 


The evolution equations for shear and Gauss curvature for the anisotropic Bianchi V space-time are given by 


a = — 


f"k 

7 r . 


2 K= --6( 2 K ) 
3 


(50) 

(51) 


The above equations are very useful in describing the evolution of shear and Gauss curvature of the 3-spheres 
with reference to time. 

Using equation (20) in the generalised Friedmann field equation (46), one can obtain 


2 _ 6f'cr 2 + 6ft (p m + p r ) + 18 p 2 f'( 2 K) + 3 (f'R — f) 

e -W- 


(52) 


In addition to the variables considered in the previous Bianchi I model, here we have considered one extra 
dimensionless variable representing the Gauss curvature, given by 


= 3 pCK) 

ee 2 


(53) 


Therefore, the variables (23) and (53) are constrained by the modified Friedmann field equation (52) as 

1 = Yj 2 T U r T — 3 K T x (54) 


Using the variables (23) and (53) along with the constraint equation (54), one can write the following 
evolution equations corresponding to the variables E, U r , K and x as: 

jyi V' 

-= —[-3 + 3E 2 + Q r + 3K - 3x- 9C(R)x + 9D(R)(E 2 + 3 K- 1)] 

dr 2 

— = + 3E 2 + n r + 3K- 3x- 9C(R)x + 9D(R)(Y? + 3 K)] 

dr 

dK 

— = K[1 + 3Z 2 + + 3K - 3x - 9C(R)x + 9D(K)(E 2 + 3 K+ 1)] 

dr 

CLX, 

— = x[3 + 3E 2 + U r + 3K - 3x + 9C(£)(1 - x) + 9D(R)(Y 2 + 3 K)] (55) 

dr 

In the dynamical system (55), it is noticed that only time (or r) dependence is contained in the parameter C 
( R ) and D(R). In the regime R f i corresponding to the Hu-Sawicki model, one can find that C(R) 1 <md D 
(R) <C 1 and hence are neglected in the analysis. Therefore the dynamical system (55) does not depend on time 
(or r) explicitly, which proves that the system we have obtained in case of Bianchi V cosmology is an 
autonomous system. 


5.1. The fixed points and exact solutions: 

The fixed points that are used to procure the exact solutions for the anisotropic Bianchi type V model are 
obtained by setting (55) equal to zero. In case of this anisotropic model we encounter anisotropic as well as 
isotropic fixed points, given by 


(E, £2 r , Kj x) = (d=l, 0, 0, 0) 

P r : (E, U r ,FC, x) = (0, 1,0, 0) 
P m : (E, U r , K, x) = (0, 0, 0, 0) 

P k : (E, O r , K, x) = jo, 

Pd : (E, a, K, x) = (0, 0, 0, 1) 
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Table 3. Fixed points and their solutions for EoS parameter, deceleration parameter, average scale factor, shear, Gauss 
curvature and physical behaviour for Bianchi V cosmology in the background of Hu-Sawicki model. 


Points 

n m 

EoS 

parameter 

( w eff) 

Deceleration 

parameter 

(q) 

Average 

Scale 

factor(n) 

Shear 

(o’) 

Gauss 

curvature 

( 2 K) 

Physical 

behaviour 

if 

0 

1 

3 

2 

a 0 \t\l 

a 0 aQ 3 \t\~ l 

0 

Decelerated 

Pr 

0 

1 

3 

1 

a 0 \t\i 

0 

0 

expansion 

Decelerated 

Pm 

i 

0 

1 

2 

a 0 \t\l 

0 

0 

expansion 

Decelerated 

Pk 

0 

1 

— 3 

0 

flobl 

0 

2 K 0 a 0 \t\ 2 

expansion 

Decelerated 

Pd 

0 

-1 

-1 

a 0 \t\> 

0 

0 

expansion 

Accelerated 








expansion 


The equation of state (EoS) parameter w e jy and the deceleration parameter q in terms of the dynamical 
variables in case of Bianchi V space-time are given by 


W eff : 


% 


% 1 - £ 2 (E 2 + 3 K) 


[e 2 + ^ 


+ 3D(R)(1 + E 2 + 3 K) - £ 2 (E 2 + K) + 2±- 

& 


+ K-( 1 + 3 C(R))x 




(56) 


q — — 1 + — 
2 


1 + H 2 + -fi r 
3 


K - (1 + 3 C(R))x + 3D(£)(X 2 


3K + 1) + 2— 

& 


(57) 


Considering C(R) <C 1 ,D(R) <C 1,^ ~ 1 and substituting the value of each variable in equations (56) and 
(57), one can obtain the EoS parameter w e jf and deceleration parameter q for the associated fixed points 
corresponding to Bianchi V model. In figure 6, we plot the evolution of EoS parameter and deceleration 
parameter along with the dimensionless dynamical variables for this anisotropic model. These figures depict that 
the Universe evolves from the anisotropic phase (w e jy = 1/3, q = 2) to the isotropic radiation dominated 
(w e ff = 1/3 ,q = 1), matter dominated (w e jy= 0, q = 1/2) and spatially non-flat (w e jy = —1/3 ,q = 0) phase, 
respectively, and enters the stable phase of late times de Sitter expansion (w e jy = — 1 ,q = —1). 

The EoS parameter w e jy and deceleration parameter q along with average scale factor a and physical 
behaviour are summarised in table 3. Table 3 shows that the points P±, P r and P m admit power law solution 
whereas P and Pj attain Milne and de Sitter solution respectively. These solutions are plotted in figure 7. 
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Figure 7. Evolution of scale factor (a) corresponding to Bianchi V cosmology in the background of Hu-Sawicki model for 
a o — = 1 - 



Figure 8. Evolution of shear (<r), Gauss curvature ( 2 K ) corresponding to Bianchi V cosmology in the background of Hu-Sawicki 
model for cr 0 = l, 2 i^o = landn 0 = 1- 


5.2. Evolution of Shear: 

Now we analyse the trace free Gauss-Codazzi equation (50), in order to find the evolution of shear for the 
anisotropic fixed points (P±). The Gauss-Codazzi equation (50), in terms of dynamical variables (23), (53) for 
Bianchi V space-time is given by 


a = 


--=[(1 + 3D(R))&]d 2 

V3 


(58) 


Integrating the above equation we obtain 


cr = a 0 a 3 = a 0 a 0 3 \t\ 1 (59) 

The above solution implies that the shear corresponding to each anisotropic fixed point diminishes with the 
evolution of cosmic time and this is displayed in figure 8. 

5.3. Evolution of Gauss curvature of the 3-spheres: 

The evolution of Gauss curvature of the 3-spheres for the spatially non-flat fixed point (i\) can be explored from 
equation (51). Integration of the equation (51) gives 
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Table 4. Eigenvalues and stabilities of the fixed points for 
Bianchi V cosmology in the background of Hu- Sawicki 
model. 


Points 

Eigenvalues [Ai, A 2 , A 3 , A 4 ] 

Stability 


[2, 3,4, 6] 

Unstable 

Pr 

[-1,1,2,4] 

Saddle 

Pm 

l-f,-1, 1,3] 

Saddle 

Pk 

[-2, -2, -1,2] 

Saddle 

Pd 

[-4, -3, -3, -2] 

Stable 


2 K = 2 K 0 a~ 2 = 2 K 0 dQ 2 \t\- 2 (60) 

The above equation shows that the Gauss curvature dies down as time progresses and this is shown in 
figure 8. 

5.4. Stability analysis of the fixed points: 

In this section, we analyze the eigenvalues of the system of equation (55), in order to study the stability 
corresponding to the Bianchi V cosmology. The eigenvalues obtained by linearising the system of equation (55) 
are as follows: 

Pf: [X u A 2 , A 3 , A 4 ] = [2 + 9D(R), 3 + 9D(R), 4 + 18 D(R), 6 + 9C(R) + 9D(R)\ 

P r : [A b A 2 , A 3 , A 4 ] = \l,4 + 9C(R), 2 + 9D(R), -1 - | D(R ) 

P m - [Ab A 2 , A 3 , A 4 ] = 1-1, 3 + 9C(R), 1 + 9 D(R), -- - —D(R) 

L 2 2 

Pm- [Ab A 2 , A 3 , A 4 ] = [-1 - 9 D(R), -2 - 9 D(R), -2 - 9 D(R), 2 + 9C(P) - 9 D(R)] 

Pd- [Ai, A 2 , A 3 , A 4 ] = [-3 - 9C(R), -4 - 9C(R), -2 - 9C(R) - 9D(R), -3 - - C(R ) - -9 D(R) 

L 2 2 

Considering C(R ) <C 1 and D(R ) <C 1, the eigenvalues regarding the Bianchi V cosmology are modified to 
the form as summarised in table 4. In this table, we also list the stability of the fixed points corresponding to the 
dynamical system. For the anisotropic fixed points P± y the linearized matrix contains four positive eigenvalues. 
Hence, these fixed points correspond to unstable nodes. The points denoting the radiation dominated, matter 
dominated and spatially non-flat universe given by P n P m and P^ respectively, are found to be saddle, as the 
linearized matrix have positive as well as negative eigenvalues. Finally, for the isotropic de Sitter fixed point all 
the eigenvalues are found to be negative, describing the stable late time accelerated expansion of the Universe. 

5.5. Phase space analysis: 

This section is devoted for the investigation of the system using phase space analysis in order to study the 
dynamical behaviour of the Bianchi V cosmology corresponding to the Hu-Sawicki f(R) model. The phase 
portrait analyses of the system for this model are displayed in figure 9. In figure 9(a) we observe the heteroclinic 
trajectories of the form 

—> Pm — > Pk Pd 

This figure is plotted for = 0, which shows the progression of the unstable anisotropic phase to the stable 
isotropic de Sitter expanding scenario following by the matter dominated and spatially non-flat isotropic age. 
Figure 9(b) shows the heteroclinic trajectories for E = 0, that gives 

P r —► Pm —> Pk — > Pd 

This figure displays the evolution of the radiation dominated phase to the matter dominated, than to 
spatially non-flat and finally to the accelerated expanding period. Figure 9(c) displays the heteroclinic 
trajectories for K = 0, that gives 



This figure depicts the progression of the unstable anisotropic universe to the stable expanding de Sitter 
phase following the radiation and matter dominated eras. 
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(a) (b) 



(c) 


Figure 9. Phase portraits for Bianchi V cosmology in the background of Hu-Sawicki model. 


6. Conclusion: 

Field equations and the fixed points for Palatini f(R) gravity in the background of anisotropic Bianchi I and V 
have been derived employing the dynamical system technique. We present the solutions associated with each 
fixed point considering an exclusive form of f(R) initially proposed by Hu and Sawicki. 

Our work begins with one of the simplest type of anisotropic model known as Bianchi type I. In this model, 
we have identified two anisotropic (P±) and three isotropic (P r , P m , P d ) fixed points. Stable de Sitter solution 
exists for the isotropic fixed point describing the late time accelerating phase of the Universe. Points indicating 
the radiation and matter dominated epochs behave as saddles whereas the anisotropic fixed points describing the 
shear are found to be unstable. Thus, in this model, the Universe starts from the anisotropic scenario that 
isotropizes with the standard radiation and matter dominated era and attain the late time de Sitter accelerating 
epoch. 

The work is further extended to another form of anisotropic model known as Bianchi type V. In addition to 
the points obtained in the previous model, this model gives rise to an additional isotropic fixed point (i\) 
describing a non-flat universe. Thus, in this case, we have two anisotropic (P±) and four isotropic (P r , P m , P^, P d ) 
fixed points. The isotropic de Sitter point describing the acceleration phase of the late universe is found to be 
stable whereas the radiation, matter and spatially non-flat isotropic fixed points behave as saddle points. The 
point corresponding to anisotropic age is found to be unstable. Therefore, this model is able to accomplish a 
sequence originating from an anisotropic phase which subsequently ends up in an isotropic de Sitter expanding 
stage, following a radiation, matter dominated as well as spatially non-flat universe. 

Anisotropic models of these types which are very effective in explaining the anisotropy of early universe, 
depict that the anisotropy dies down as (time) -1 . The explicit treatment of the anisotropic model Bianchi type V 
furnishes the Gauss curvature evolution and shows that curvature dissipates as (time) -2 . 

In [5], the authors have been able to obtain the sequence of radiation-dominated, matter-dominated and late 
time acceleration for FLRW cosmology considering f(R ) = R — /3/R n . Authors of [11] could reproduce the 
same sequence considering the Hu-Sawicki/(R) gravity in Palatini formalism. Working in the background of 
anisotropic Bianchi I and V model, we have been able to regenerate the same sequence. Additionally, we have 
been able to obtain a sequence that starts from anisotropic universe which ends into a stable accelerating phase. 
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The authors of [28-30, 39, 44-47] have shown that the anisotropic universe accomplished the isotropic stage 
with the evolution of time in case of anisotropic Bianchi models, which is consistent with our analysis. 
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